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ABSTRACT 

To date, several exoplanets have been discovered orbiting stars with close binary companions (a < 30 
AU). The fact that planets can form in these dynamically challenging environments implies that planet 
formation must be a robust process. The initial protoplanetary disks in these systems from which 
planets must form should be tidally truncated to radii of a few AU, which indicates that the efficiency 
of planet formation must be high. Here, we examine the truncation of circumstellar protoplanetary 
disks in close binary systems, studying how the likelihood of planet formation is affected over a range 
of disk parameters. If the semimajor axis of the binary is too small or its eccentricity is too high, 
the disk will have too little mass for planet formation to occur. However, we find that the stars 
in the binary systems known to have planets should have once hosted circumstellar disks that were 
capable of supporting planet formation despite their truncation. We present a way to characterize 
the feasibility of planet formation based on binary orbital parameters such as stellar mass, companion 
mass, eccentricity and semi-major axis. Using this measure, we can quantify the robustness of planet 
formation in close binaries and better understand the overall efficiency of planet formation in general. 
Subject headings: planets: formation — protoplanetary disks — binaries: close 


1. INTRODUCTION 

Close binary systems are a challenging environment 
for planet formation because a binary companion will 
dynamically perturb the protoplanetary disk. However, 
several exoplanets have been found in S-type orbits (cir¬ 
cumstellar) around stars in binary sys tems with semi- 
majo r axes of less than 30 AU (e.g. iRaghavan et al.l 
1200611 . These systems represent interesting test cases 
for planet formation theories. By comparing the fre¬ 
quency of planets in binary systems with that around 
single stars, we can learn about how robust the planet 
formation process can be. In this paper, we address the 
topic of circumstellar planets (S-type orbits), as opposed 
to circumbinary planets (P-type orbits). 

HD 188753 is a hierarchical triple system, with a 12.3 
AU separation and orbit al eccentricty o f 0.5 between the 
A and BC components. iKonackil (|2005[ 1 claimed the de¬ 
tection of planet around the A component. However, we 
found that the binary parameters for HD 1 88753A do not 
allow for in situ planet formation to occur (iJang-Condelll 
120071 henceforth Paper I). Indeed, it was later found that 
there was no planet there after all (lEggenberger et ^ 

|2003). 

On the other hand, 7 Cep has a separation of 20 
AU and eccentricity of 0.4, and has a well-established 
gas giant planet or b iting the A comp onent at ~2 AU 
(jHatzes et al.l 120031 : lEndl et ^ 1201111 . Although at 
first glance, its orbit is not very different from that of 
HD 188753, we found that the binary parameters of 

J ' Cep A do allow in si tu planet formation to occur 
Jang-Condell et aLll2008L henceforth Paper 11). 

Circumstellar planets have been discovered in a few 
other binary systems of comparable sepa r ations. These 
include HD 196885 (iCorreia et al l 120081 iFischer et al.l 
120091 : IChauvin et al.ll2011[ l. a Cen (|Pourbaix et al.ll2002l : 

hjangcon@uwyo.edu 


iDumusoue et al.lF2012HHatzesll2013[ l. and HD 41004. In¬ 
terestingly, HD 41004 has substellar companions in both 
A and B componen ts: HD 41004 Ab is a 2.5 Mj planet 


(Zucker et al. 


Zucker et al 


2004) and HD 41004 Bb is a brown dwarf 


200311 . The binary orbits and planet char¬ 


acteristics of the systems discussed here are listed in Ta¬ 
ble [TJ 

One way to study the robustness of planet formation in 
binaries is to study the long-term dyn amical stability of 
those planets iHolman fc WiegertI le.g. Il999l l. However, 
the focus here in this work is on the formation of planets 
in the first place. Rather than doing a detailed calcula¬ 
tion ofthe_^namics of planet formation as in, for exam¬ 
ple, [Thebaula (|2011l l. we want to understand planet for¬ 
mation more broadly by considering the effects of binary 
interactions with the initial protoplanetary disk. That is, 
given that the protoplanetary disk will be truncated by 
a binary companion, what is the likelihood that planet 
formation can take place at all? 

Calculations of disk truncation by binary companions 
show that stars in b inary systems are still able to retain 
circumstellar disks (lArtvmowicz fc Lubowl [l994f) . This 
is supported by millimeter and submillimeter observa¬ 
tions of pre-main sequence binaries which show that bi¬ 
nary stars with separations of tens of AU are able to re¬ 
tain their disks iTensen et'ffilllQQfiHOsterloh fc BeckwithI 


119951 : [Jensen et al.lll994 iBeckwith et al.lll990ll . 

The objective of this paper is to give a measure for the 
difficulty of planet formation across a range of binary 
star parameters. We consider S-type orbits, in which the 
planet orbits one of the stars in the binary system, as op¬ 
posed to P-type orbits, which refer to circumbinary plan¬ 
ets. We examine the initial protoplanetary disk around 
a given star and calculate the degree of truncation of 
this disk. Rather than determining where exactly plan¬ 
ets should form in these truncated disks, we focus on 
whether planet formation is possible at all. 
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TABLE 1 

Binary systems analyzed in this paper. 


HD 188753A 

7Cep A 

HD 41004A 

HD 41004B 

HD 196885A 

«Cen B 

M* (Mq) 

1.06 

1.40 

0.7 

0.4 

1.31 

0.93 

Mb (Mq) 

1.63 

0.41 

0.4 

0.7 

0.45 

1.1 

n 

0.39 

0.77 

0.64 

0.36 

0.74 

0.46 

a (AU) 

12.3 

20.2 

20 

20 

21 

23.5 

e 

0.5 

0.41 

0.4 

0.4 

0.42 

0.52 

nip sini (Mj) 

— 

1.85 

2.5 

18.4 

2.96 

0.0034 

Up (AU) 

— 

2.1 

0.006 

7 X lO-'^ 

2.6 

0.04 

ep 

— 

0.05 

0.7 

0.08 

0.46 

0 

reference 

a 

b 

b,c 

d 

e 

f 

Nm 

9 

12 

12 

12 

12 

12 

^D1 

0 

1 

1 

1 

1 

1 

^CA 

0 

8 

7 

4 

8 

3 


Note. — Orbital parameters without subscripts (a and e) refer to those of the binary. The subscript ‘p’ refers to planetary or substellar 
companions. 

Ref erences . — falTE ggenberger et ^ Il200ff) . ('bi lEndl et al.l 11201111 . Icl IZucker et al.l II2004I ). Idl lZucker et al.1 11200311 . lel lCorreia et aO 
II2008I ). and ff) IDumusque er'an ~ ll^l2l 1 


In this paper, we extend the methods of determining 
the feasibility of planet formation in close binaries as de¬ 
veloped in Papers I and II to a wider range of binary 
star configurations. For a given set of binary parame¬ 
ters, we develop a formalism for describing the feasibility 
of planet formation as measured by the number of disk 
models that would allow for either core accretion or disk 
instability to occur. 

2. METHODS 

2.1. Truncated Disk Models 

We base our estimates of the likelihood of planet forma¬ 
tion for a given binary system on the range of disk param¬ 
eters that permit enough mass in the disk for planet for¬ 
mation to occur, given that a circumstellar protoplane¬ 
tary disk will be truncated by a binary companion. These 
methods are explained fully in Papers I and II, and are 
summarized briefly here. 

For convenience, we refer to the planet-hosting star 
as the primary and the binary stellar companion as the 
secondary, regardless of the relative masses between the 
binary components. Accordingly, M* and Mb are the 
masses of the primary and secondary, respectively. The 
mass ratio of the binary is defined as /r = M*/ {M^ + Mb). 
The binary orbit is defined by the semimajor axis a, and 
eccentricity e. 

The primary mass, M*, affects both the gravitational 
potential of the disk as well as the amount of stellar irra¬ 
diation of the disk since luminosity varies strongly with 
stellar mass. We consider a range of masses for the pri¬ 
mary star from 0.5 to 2.0 Mq. To calculate the lumi¬ 
nosity of the primary star, needed as an input parameter 
for the disk models, we use pre-main seq uence models for 
stars of I Myr of age (jSiess et al.ll2000l i. The values for 
Teff, R* and L, for each value of M* selected are listed 
in Table [21 

The two additional parameters that determine the 
properties of a disk are the mass accretion rate (M) and 
viscosity parameter (a). For each stellar mass, we calcu¬ 
late a suite of disk models by varying a and M, as 

a G {0.001,0.01,0.1} (1) 

and 

M G (10-®, 10"®, 10"^ 10"®, 10"®, 10-^1 M0 yr-^ (2) 


TABLE 2 

Stellar parameters as a function of stellar mass 

Stellar mass Tgg R* LD 
(Mr) (K) (Rr) (La) 


0.5 

3770 

2.1 

0.8 

0.7 

4026 

2.4 

1.4 

0.9 

4209 

2.5 

1.8 

1.0 

4282 

2.6 

2.0 

1.1 

4343 

2.7 

2.3 

1.3 

4447 

2.9 

3.0 

1.5 

4536 

3.1 

3.7 

1.7 

4616 

3.3 

4.4 

1.9 

4690 

3.4 

5.0 

2.0 

4721 

3.5 

5.5 


fDerived from and i?*. 


for a total of 18 disk models for each set of binary pa¬ 
rameters a, e). The lower accretion rates reflect 

values typically observed in T Tauri disks. The highest 
accretion rates are more typical of active FU Orionis-like 
objects or very young protoplanetary disks than of pas¬ 
sive T Tauri disks. Nevertheless, we include these high 
values of accretion for completeness. 

A given value of a sets the viscosity, according to 

V = aCsH (3) 

where Cs is the ther mal sound speed and H is the scale 
height of the disk (iShakura fc Sunvae^ 119731) . These 
quantities are related hj H = CgD. where fl is the Ke- 
plerian angular speed, and Cs = ^/kT/p. where k is the 
Boltzmann constant, T is the midplane temperature, and 
p is the mean molecular weight of the gas, whose assumed 
composition is that of molecular hydrogen. The viscosity 
sets the rate of the mass flow through an accretion disk, 
according to 

izF = ^(l-r-/A.) (4) 

where E is the surface density, r is the s tellocentric ra - 
dius, and i?* is the stellar radius (see, e.g.. lPrin ilifTMl) . 
If we assume that the accretion rate, M, is constant with 
r, then Eq.Ujgives a relation between T and E in the disk. 
Assuming that the only heating sources in the disk are 
viscous accretion and stellar illumination, we calculate 
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Fig. 1.— Disk profiles for the least massive (0.5 Mq, left) 
and most massive (2.0 Mq, right) planet host stars examined 
in this work. We show the total disk mass (top), total dust 
mass (middle), and Toomre Q parameter (bottom) versus radius 
for the disk parameters explored. The accretion rate (M) is 
indicated by line color, while values of the viscosity parameter of 
a = 0.001, 0.01, and 0.1 are indicated by solid, dotted, and dashed 
lines. The truncation radius for each disk for an equal-mass binary 
companion = 0.5) at a = 30 AU is shown by the vertical gray 
line for e = 0, by hollow symbols for e = 0.3 and filled symbols for 
e = 0.6. 


T(r) and E(r) self-consistently. Thus, a and M can be 
selected to “tune” the mass of the disk. 

The truncation radius for each of the 18 disk models 
for a selected primary mass is dete r mined from values 
calculated in lArtvmowicz fc Lubowl (jl994[ i , who define 
the truncation radius to be where viscous torques bal¬ 
ance the tidal torques. This truncation radius depends 
both on the Reynolds number, which is controlled by the 
viscosity v, and the binary orbit, represented by the mass 
ratio and eccentricity. That is, for each set of M and a 
for the disk, and M*, /i, a, and e for the binary, there is 
a unique disk truncation radius. 

In Figure[TJ we show how the structure of the disks vary 
with input parameters. The two sets of plots show disk 
profiles for the least massive (0.5 Mq) and most massive 
(2.0 Mq) stars explored. Each line represents a selected 
value of M and a for a given disk model, as indicated by 
color and line type according to the legend. For each disk 
model, we show the enclosed disk mass (top), total dust 
mass (middle), and Toomre Q parameter (bottom) as a 
function of radius. The Toomre Q parameter is given by 


Q 


CsK 

ttGS’ 


( 5 ) 


where k is the epicyclic frequency and G is the gravita¬ 
tional constant. Since the disk orbital velocities are close 
to Keplerian, we assume that k Ki fl. 

If the binary orbit is circular, then all disks will be 
truncated at the same radius. The truncation radius for 
an equal-mass binary companion (/r = 0.5) at 30 AU in 


a circular orbit is 10 AU, as indicated by a gray line in 
Figure [T] However, if the binary companion is eccentric, 
the truncation radius varies with disk properties. This 
is because the Reynolds number is different in each disk, 
therefore the balance of tidal and viscous torques occurs 
at different radii. The truncation radii for a binary orbit 
with /i = 0.5 and a = 30 AU with e = 0.3 (0.6) is 
shown by hollow (filled) symbols in Figure [1] As the 
eccentricity increases, the truncation radii move inward 
as expected, since the periastron occurs at smaller radius 
with increasing eccentricity. 

The disk profiles shown in Figure [T] vary more with M 
and a than with stellar mass. As the stellar mass in¬ 
creases, the disks become slightly more massive overall, 
but with less dust and higher values of Q. This is be¬ 
cause the increased stellar illumination heats the disks to 
higher temperature, which affects both disk sublimation 
and pressure support of the gas. 

The Toomre Q parameter decreases with radius in 
most models of protoplanetary disks as a general rule. 
In the disk models presented, the midplane temperature 
and suface density typically behave as oc and 

E oc respectively. Since H oc the overall 

variation of Q is to decrease with r, 

2 .2. Core Accretion vs. Disk Instability 

Given a set of parameters defining a particular binary 
system (M*, /i, a, and e), we can count how many of the 
18 truncated disk models constructed around the pri¬ 
mary star ( 1 ) have sufficient total mass to form planets, 
( 2 ) have sufficient solid material to form a planet core, 
or (3) have low enough value of Q for disk instability to 
occur. If either condition (2) or (3) is satisfied, then (1) 
is also satisfied, so we henceforth consider only the more 
stringent criteria for planet formation based on core ac¬ 
cretion and disk instability, rather than the total disk 
mass. We define the number of disks in which core ac¬ 
cretion is possible based on the mass of solid material to 
be Nca and the number of disks in which disk instability 
is possible based on the value of Q to be A^di- 

In the core accretion process of planet formation, gi¬ 
ant planets are formed by first creating a rocky core from 
the agglomeration of dust particles in the disk. When it 
reaches sufficient mass, typically around 10 Mq, it be¬ 
comes massive enough to accrete and retain a gaseous 
atmosphere, and gas giant planet is born. Our criterion 
for determining that a given truncated disk model is ca¬ 
pable of forming giant planets by core accretion, then, is 
that it contains at least 10 Mq of solid materials, i.e. dust 
and ice. If the total amount of dust in a given truncated 
disk model is above 10 Mq, then we say that that disk 
is capable of supporting planet formation by core ac¬ 
cretion. So, Nca is the total number of truncated disk 
models that can support planet formation by core accre¬ 
tion for a given set of binary parameters. For example, 
Nca for M, = 0.5 Mq, /i = 0.5, a = 30 AU, and e = 0 
can be read off the dust mass plot in Figure [T] (left cen¬ 
ter) by counting the number of profiles that intersect the 
truncation radius (gray line) at 10 Mq or above, giving 
Nca = 12. When the eccentricity is raised to e = 0.3, the 
solid content of the disk model with M = 10“®MQyr“^ 
and a = 0.1 drops below 10 Mq, so Nca = H- 

The criterion for planet formation by disk instability, 




















4 


JANG-CONDELL 


on the other hand, is determined by the value of the 
Toomre instability parameter, Q. When Q < 1, the disk 
is gravitationally unstable and can rapidly form gas giant 
planets through fragmentation. If a given truncated disk 
model has Q < 1 at any point in the disk, then we de¬ 
clare that the disk is capable of forming planets by disk 
instability, and A^di is the total number of truncated disk 
models that are gravitationally unstable for a given set 
of binary parameters. For M* = 0.5 Mq, /r = 0.5, a = 30 
AU, and e = 0, AA)i = 4, and when the eccentricity is 
raised to e = 0.3 IVdi = 3. 

3. RESULTS AND ANALYSIS 

In Papers I and II, we used a fixed set of parameters for 
M*, a, e, and ^ and used those parameters to calculate 
Nca, and A^di- Here, we extend the ranges of A/*, a, e, 
and /i in order to see how Nqa and A^di vary with each 
parameter. In TablelU we show our calculations for Nca 
and A^di for the binary systems listed there. 

Our previous results for HD 188753A (Paper I) and 7 
Cep A (Paper H) are shown in the first two columns, re¬ 
spectively. We also repeat our calculations for the close 
binary systems HD 41004, HD 196885A, and a Cen B. 
Although HD 41004Bb is a brown dwarf rather than a 
planet, we include that system in this paper as a com¬ 
parison. 

HD 188753A, where there is no confirmed planet, is the 
only system with A^di = 0 ; all other systems have A^di = 
1. Similarly, HD 188753A is the only system examined 
here with Nqa = 0 . If we use Nqa as a proxy for the 
difficulty of planet formation, then a Cen B should be 
the least likely system to harbor a planet, and 7 Cep A 
and HD 196885A the most likely. 

In Figure [5] we show examples of how Nqa and A^di 
vary with a and e, holding M* and fi fixed. That is, the 
masses of the binary system are held fixed and their or¬ 
bits are allowed to vary. One plot shows where 7 Cep A 
falls in parameter space, while the other shows the loca¬ 
tion of HD 188753A. As eccentricity decreases and semi¬ 
major axis increases, Nqa (green solid contours) and A^di 
( red dashed contours) both increase, as expected. This 
is because as periastron increases, the disk truncation 
radius also increases, allowing a larger disk to remain 
around the star. A^di increases also because Q tends to 
decrease with increasing distance from the star, and Nqa 
increases because more cold disk material remains in the 
disk. 

The irregularity in the spacing of the green contours 
(Nca), particularly at large a and small e, results from 
the finite sampling of disk models. If we had included 
more values of M and a, then there would be more con¬ 
tours in this area. Some of the contours bunch up at 
small separations and low eccentricities, largely because 
the solid content of the disks are not a simple function 
of M and a, as shown in the middle panels of Figure [TJ 
The dust mass versus radius profiles of the disks cross 
each other at small separations, with steeper profiles at 
higher accretion rates. This is because holding a con¬ 
stant, higher accretion rates give higher disk mass ac¬ 
cording to Equation but the heating from this accre¬ 
tion causes sublimation of the refractory species. These 
competing effects mean that at large disk radii, the total 
solid content should overall be higher in the more mas¬ 
sive, higher accretion rate disks, but at small radii, the 


refractory species sublimate and those disks have lower 
solid content. 

The values of A^di and Nqa can be fit to a polynomial 
function dependent on a, e, /r, and M*. The general form 
of this polynomial is 


N — ^ ^ ^ijkl 


Vl AU/ 





( 6 ) 


Since we want to keep the analytic function as simple 
as possible, we do not try to fit the bunched contours, 
but rather try to find a smooth fit that will give a rea¬ 
sonable approximation for the true values of Nca and 
A^di- We find best-fitting polynomial coefficients, Cijki, 
for A^ca, using those points where they have tabulated 
values where 0 < Nca < 10 , because this avoids the 
plateau in Nca values at larger separations and low 
eccentricity. For A^di, we use tabulated values where 
0 < A^di < 5. As shown in Figure [U values of A^di higher 
than 5 occur only at a > 50 AU and at low eccentricity, 
and in a region where Nca > 10 . 

We find that for Nca we get reasonably good fits for 
a function linear in a and e and quadratic in and M*. 
For Ndi, the fitting function is linear in all four parame¬ 
ters. The best-fitting coefficients are tabulated in Table 
[21 Comparing the directly calculated and tabulated val¬ 
ues for Nca to the polynomial fit for all tabulated values 
where Nca < 10 , the root-mean-square deviation is 0.62. 
For A^di < 5, the rms deviation is 0.44. 

The values of Nca and Nui for the binary systems 
considered in this paper as calculated using this polyno¬ 
mial equation are tabulated in Table 01 Since both Nca 
and A^di cannot be less than 0 , negative values should 
be interpreted as 0. The resulting values are all within 
±1 for Nca, and within ± 0.6 for A^di of the directly cal¬ 
culated values shown in Table 01 The polynomial values 
are more discrepant for larger values of Nca, underesti¬ 
mating the true values. This is a result of the polynomial 
equation smoothing over the bunching of Nca contours, 
as discussed above. 

The values of A^di for all binary systems considered 
are less than 2, as shown both in Table 01 and Table 01 
This could be an indication that core accretion is more 
likely than disk instability in close binaries. As shown 
in Figure 0 ] the disk model with the lowest Q values has 
the highest accretion rate, 10 “'^ Mq yr“^ and the lowest 
viscosity parameter explored. Disks with this high an 
accretion rate are typically outburst systems, such as FU 
Orionis objects. Accretion rates this high are typically 
transient. The disk contains 100 Mj or 0.1 Mq within 
just 1 AU, so it is extremely massive, hence it is the 
most gravitationally unstable. Because a steady-state 
disk with these properties is unlikely, planet formation 
via disk instability is also unlikely. 

If we were to set a minimum threshold for Nca for 
which planet formation can occur in a close binary given 
the binary systems explored in this paper, we first need 
to consider which systems to include. The existence of a 
Cen Bb has yet to be definitively confirmed, and it has 
the lowest value of Nca of 3, not including HD 188753A. 
If a Cen Bb is confirmed, then we might set A^cA.mm = 3. 
One might dispute that HD 41004Bb is a brown dwarf 
and ought not be considered, also. Then, we might 
set A^cA.min = 7, since 7 Cep A, HD 41004A, and HD 
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semi-major axis (AU) semi-major axis (AU) 


Fig. 2.— Planet formation likelihood plots for HD 188753A (left) and 7 Cep A (right). Each plot shows how A’qa -^DI vary with 
eccentricity and semi-major axis for the binary orbit for a fixed primary and secondary mass. Solid magenta contours show values for A’qa 
and dashed blue contours show those for A’di- The positions of HD 188753A (left) and 7 Cep A (right) in a-e space are indicated by 
crosses. 


TABLE 3 

Coefficients for Equation (O. 




II 

p 

II 

O 

i = l,j = 0 

i = 0 , j = 1 

i = l,j = l 

Core Accretion 

II 

o 

1 = 0 

-6.7585 

0.4308 

2.9369 

-0.4411 

k = 1, 

1 = 0 

-2.0956 

1.6243 

-9.5050 

-1.4578 

k = 2, 

1 = 0 

1.8452 

-0.5589 

16.0528 

0.3879 

II 

o 

1 = 1 

1.5486 

-0.1609 

-1.9967 

0.1696 

k = 1, 

1 = 1 

-9.3964 

0.2216 

3.5221 

-0.3374 

k = 2, 

1 = 1 

7.9117 

-0.3244 

-5.4470 

0.4971 

II 

o 

1 = 2 

-0.2951 

0.0364 

0.4372 

-0.0363 

k = 

1 = 2 

2.7651 

-0.0993 

-1.6003 

0.1275 

k = 2, 

1 = 2 

-2.2214 

0.0920 

2.1346 

-0.1406 

Disk Instability 

fc = 0 , 

1 = 0 

0.3921 

0.0521 

-0.8329 

-0.0431 

fc = 1 , 

1 = 0 

-0.1135 

0.1106 

1.1926 

-0.0975 

o' 

II 

1 = 1 

-0.2704 

-0.0066 

0.3131 

0.0012 

k = 

1 = 1 

0.2944 

-0.0296 

-0.7153 

0.0327 


TABLE 4 

Calculated values of and Nj^i from Equation (O. 



HD 188753 

yCep 

HD 41004A 

HD41004B 

HD 196885 

CfCenB 

Ndi 

0.54 

1.41 

1.57 

1.25 

1.45 

1.15 

No A 

-2.19 

7.04 

6.87 

3.83 

7.13 

3.45 


196885A all satisfy this criterion. Setting VcA.min = 7 is 
the more stringent constraint for planet formation, since 
we require more disks to satisfy the requirement for suffi¬ 
cient solid mass. Based on these considerations, we pre¬ 
dict that Nca ^ 7 for planet formation to be readily 
feasible, based on the known planets in binary systems. 
For 3 < Nca < 7, planet formation may still be possible, 
as in a Cen A, but less common. 

In Figure [3J we show the allowed parameter space for 
planet formation in close binaries, with VcA,min = 3 
(left) or VcA.min = 7 (right). That is, we set equation 
([6|) equal to 3 or 7, choose fixed values of M* and and 
draw the resulting curve in a-e space. Since equation m 
is linear in a and e, the resulting curve for N = VcA.min 


is hyperbolic in a-e space for fixed M* and /i. 

For a given binary pair, giant planet formation by core 
accretion is allowed for greater separations and lower 
eccentricity than the indicated contour, or under each 
curve. Planet formation also becomes more feasible as 
the primary mass decreases and the mass ratio increases, 
as shown by the shifting of the contours up and to the 
left. 

Note that for a = 12 AU, the semi-major axis of the 
HD 188753 binary, planet formation is highly disfavored. 
Even considering A'cA.min = 3, only stars with relatively 
low-mass companions {fi > 0.7 or M^/M* < 3/7) and 
low eccentricity (e < 0.5) might allow planet formation 
to occur. However, for a = 20 AU, a large range of binary 
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Fig. 3.— Contours in a-e space for = 3 (left) and AIqa = 7 (right) for varying binary masses. Each contour represents a selected 

choice for M, and fi, and the axes show allowed values of a and e. The mass of the primary star is labeled by color, and the mass ratio 
of its binary companion is indicated by line style: dot-dashed, dashed, dotted, and solid for fi = 0.1, 0.3,0.7, and 0.9, respectively. Planet 
formation becomes more feasible to the right and below each curve. 


properties allow for planet formation to occur, although 
larger values of fj, are favored. 

4. DISCUSSION AND CONCLUSIONS 

We have calculated the feasibility of giant planet for¬ 
mation in the close binary systems HD 188753, 7 Cep, 
HD 41004, HD 196885, and a Gen. We find that ex¬ 
cept for HD 188753, where there is not a planet, core 
accretion is the more likely planet formation mechanism 
than disk instability based on simple models of the trun¬ 
cated protoplanetary disk from which the planet formed. 
These truncated disks contain a sufficiency of solid ma¬ 
terial from which planet cores can form, but do not meet 
the criterion for gravitational instability. Based on the 
results for these systems, we conclude that for binary 
systems of separations <20 AU and eccentricity >0.4, 
circumstellar giant planets may only form via core accre¬ 
tion. 

More generally, we extend our results to binary sys¬ 
tems of arbitrary masses, separations, and eccentricities 
to predict the likelihood of giant planet formation by ei¬ 
ther core accretion or disk instability. These likelihoods 
are represented by the quantities A’ca and iVni, respec¬ 
tively. We find an empirical fit to A^ca and A^di as a 
function of binary mass, mass ratio, separation. By com¬ 
paring A^ca and A^di as predicted in this paper to the 
observed frequency of planets in binaries, we can gain 
insight into the planet formation process in general. 

The estimates of planet formation likelihood repre¬ 
sented by A^ca and Nqi are not meant to be definitive 
relations for the probability of planet formation. Rather, 
these values should be used to represent the relative ease 
with which planet formation can take place in binary 
systems. As more planets are found in close binary sys¬ 
tems, those statistics can be compared to the likelihoods 
calculated in this work to yield a better understanding 
of how planets can form in the dynamically challenging 
environment of a close binary star. 

We have considered only giant planet formation. In 
principle, the requirement for terrestrial planet forma¬ 
tion is less stringent, since less solid material is required. 
These calculations would be a possible direction for fu¬ 
ture work. 

Our results have been calculated under the assumption 


that the circumstellar disk is coplanar with the binary or¬ 
bit. It is conceivable that the disk could be misaligned 
with the binary, as is t he case for the wide pr e-main 
sequence binary HK Tan (iJensen fc Akesonll20I4D . How¬ 
ever, tidal effects from close stellar companions with < 20 
AU separations will probably cause the disk to become 
aligned. The orbits of circumbinary planets discovered 
by the Kepler mission are closely aligned with those of 
the binary, suggesti ng that disks are l ikely to be aligned 
with binary orbitsjlDovle et al.l[ 20 Tlt iWelsh et al.ll 2012 l : 

lOrosz et al.ll20I2bll^ ^ 

We have made the assumption in the core accretion 
case that likelihood of planet formation depends only on 
the total solid mass being above a certain threshold, and 
that it does not depend on the amount above that thresh¬ 
old. We do not calculate whether or not these solids can 
come together to form a single core, nor do we consider 
the possibility of the formation of multiple cores. These 
types of calculation would involve detailed modeling of 
the dynamics of the particles within the disk, and is be¬ 
yond the scope of this paper. 

We also do not consider the growth and evolution of gi¬ 
ant planets within the disk beyond the initial formation. 
Dynamical effects such as tidal torques on planetesimals 
and planet embryos, planet-planet scattering, and disk 
migration will affect the survivability of planets within 
the system. These effects have be en considered i n some 
individual systems by others fe.g. lThebaultl[2Mll) . 

We have not considered how metallicity or chemical 
composition of the disk might affect our results. In¬ 
creased metallicity will affect both the disk opacities and 
the amount of solids in the disks. Qualitatively speak¬ 
ing, increased opacity will result in slightly higher disk 
temperatures, which will decrease the refractory content 
in the most severely truncated disks. On the other hand, 
higher metallicity will raise the refractory content overall, 
so these effects will compete with each other. This could 
result in a lower likelihood of giant planet formation in 
the most tightly bound binary systems, but increased 
likelihood in wider systems. 

Thanks go to Eric Jensen and an anonymous referee for 
helpful comments that greatly improved this paper. The 
author acknowledges support from NASA ATP Grant 
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